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ABSTRACT 

We report on an analysis of RXTE data of the transient neutron star low-mass X-ray binary (NS- 
LMXB) XTE J1701-462, obtained during its 2006-2007 outburst. The X-ray properties of the source 
changed between those of various types of NS-LMXB subclasses. At high luminosities the source 
switched between two types of Z source behavior and at low luminosities we observed a transition 
from Z source to atoll source behavior. These transitions between subclasses primarily manifest 
themselves as changes in the shapes of the tracks in X-ray color-color and hardness-intensity diagrams, 
but they are accompanied by changes in the kHz quasi-periodic oscillations, broad-band variability, 
burst behavior, and/or X-ray spectra. We find that the low-energy X-ray flux is a good parameter to 
track the gradual evolution of the tracks in color-color and hardness-intensity diagrams, allowing us to 
resolve the evolution of the source in greater detail than before and relate the observed properties to 
other NS-LMXBs. We further find that during the transition from Z to atoll, characteristic behavior 
known as the atoll upper banana can equivalently be described as the final stage of a weakening Z 
source flaring branch, thereby blurring the line between the two subclasses. Our findings strongly 
suggest that the wide variety in behavior observed in NS-LXMBs with different luminosities can be 
linked through changes in a single variable parameter, namely the mass accretion rate, without the 
need for additional differences in the neutron star parameters or viewing angle. We briefly discuss the 
implications of our findings for the spectral changes observed in NS LMXBs and suggest that, contrary 
to what is often assumed, the position along the color-color tracks of Z sources is not determined by 
the instantaneous mass accretion rate. 

Subject headings: accretion, accretion disks — stars: individual (XTE J1701-462) — stars: neutron 
— X-rays: stars, binaries 



1. INTRODUCTION 

Neutron star low-mass X-ray binaries (NS-LMXBs) are 
systems in which a low-magnetic-field neutron star ac- 
cretes matter from a low-mass companion star through 
Roche-lobe overflow. They display a wide range of X-ray 
spectral and variability properties, both as a class and 



whi 



shown in Figure IT] The Z source tracks typically show 



as individual sources (see van der Klis|[2006 for a recent 
review) . Based on their correlated X-ray spectral and 
rapid variability behavior, NS-LMXBs are commonly di- 
vided into two subclasses, the so-called Z sources and the 



atoll sources ( Hasinger fc van der Klis|1989 1, named after 
the shapes they trace out in X-ray color-color (CD) and 
hardness- intensity diagrams (HID). The Z sources have 
luminosities close to the Eddington luminosity (~0.5- 
1 LEdd or more), whereas atoll sources have luminosities 
in the range ~0. 01-0. 5 LEdd- 

An example of the variety in shapes that is observed 
in the CD/HID tracks of the Z and atoll subclasses is 
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three branches, which from top to bottom are called the 
horizontal branch, the normal branch, and the flaring 
branch. Based on the presence and orientation of these 
branches, the Z sources can be further divided into the 
'Cyg-like' Z sources (Cyg X-2, GX 5-1, and GX 340-^0) 
and the 'Sco-like' Z sources (Sco X-1, GX 17-^2, and GX 
349+2), with the former showing 'Z'-shaped tracks in the 
HID (Fig.flk) and the latter more V-shaped tracks (Fig. 
[^,c; see also [Kuulkers et al.||1994| |Homan et "aI1|2007b 
(hereafter H07)). The atoll sources can show 'Z'-shaped 
tracks as well (Fig.[I^), although the individual branches 
are thought to have a different physical nature than the 
Z source branches (Barret & Olive 2002; van Straaten 
et al.||2003| |Reig et al. 2004^ van der Khs 2006) . T'he 
branches of the Z^shaped atoll tracks are called, from top 
1 



to bottom (Fig. [1^), the extreme island state, the island 
state, and the lower and upper banana branches; these 
branches are also frequently referred to as the hard, tran- 
sitional/intermediate, and soft states, respectively. The 
low-luminosity atolls are typically found in the (extreme) 
island state, whereas the high-luminosity ones tend to be 
found mainly on the (upper) banana branch (Fig. IIH). 

In addition to motion along atoll or Z branches, some 
sources also show changes in the shape and location of 
their tracks in the CD and HID, which is often referred 
to as secular motion. Examples of this can be seen in 
Figure [T] with the horizontal branch of GX 17-1-2 and 
the upper-banana branches of GX 9-1-9 and 4U 1735- 
44 being traced out at varying intensity levels. Much 
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Figure 1. Example CDs (left) and HIDs (right) of the Z source 
NS-LMXBs GX 5-1, GX 17+2, and GX 349+2, and the atoll 
source NS-LMXBs GX 9+9, 4U 1735-44, 4U 1636-53, and 4U 
1705-44. Each data point represents a 256s average. The dia- 
grams are taken from a compilation by Fridriksson et al. (2010, in 
prep.) and were created following the steps described and referred 
to in ij2]and jjs] Z and atoll branches/states are indicated in the 
CD or~HID oi each source. They are: horizontal branch (HB), 
normal branch (NB), flaring branch (FB), upper-banana branch 
(UB), lower-banana branch (LB), island state (IS), and extreme 
island state (EIS). For the Z source GX 5-1 we have also labelled 
the upturn from the horizontal branch (HB-upturn), and the 'dip- 
ping' flaring branch (dipping FB). 



stro nger secular motion ca.n be found in the Z sou rce Cyg 
X-2 ( [Kuulkers et~aI]|1996[[Wijnands et al.||1997[) , where 
it occasionally results m switches between ( Jyg-like and 



Sco-li ke Z source behavior (see Figure 1 in Muno et al. 
2002 1 . Finally, we note that transitions between the atoll 



(extreme) island state and the banana branches, and vice 
versa, often show strong hysteresis in HIDs, resulting in 
jumps between those branches (e.g., Maitra fc Bailyn 



2004|JGladstone et al.|[2007||Bellon et al ||!^0()7[ see also 
Fig. U',g). It is not clear if and how this behavior is 
related to the secular motion. 

In addition to differences in X-ray luniinositiy and 
CD/HID tracks, the NS-LMXBs subclasses also show 
differences in X-ray spectral and va riability properties 
( van der Klis| 2006 1 , radio behavior ( Migliari & Fender 



200t)[), as well as m the rate at which they show typ e 



1 (thermonuclear) X-ray bursts (Galloway et al. 2008J) 



This disparity in the properties between the subclasses 
was a long-standing problem in the understanding of the 
accretion processes in NS-LMXBs. While differences in 
the mass accretion rate (M) could explain some of the 
variety in observed properties (e.g., in X-ray luminos- 
ity and bursting behavior), it was unclear if one could 
link the wide variety in CD/HID tracks, rapid variabil- 
ity, and broadband spectral properties through changes 
in M alone. Additional parameters have been proposed 
to explain some of the observed differences, not only be- 
tween Z and atoll subclasses, but also between the Cyg- 
like and Sco-like Z sources: magiietic field strength and 
viewin g angle ( e.g., Psaltis et al. 19951 'Hasinger & van 



der Klis 1989; Kuulkers & van der Klis 1995) 



I'esting for the necessity of such additional parameters 
is complicated by the fact that none of the classical Z 
sources has been observed to venture into the atoll lu- 
minosity range (and vice versa). Until recently only two 
sources were suspected to have made a transition be- 
tween Z-like and atoll-like behavior: Cir X-1 and XTE 
J1806-246 (= 2S 1803-245). Ch X-1 was found to show 
Z-like properties at high luminosities (Shircy et al. 19991 
and atoll-lik e properties at low luminosities (,UosterbroeK| 



et al.||1995[ ). XTE J1806-246 showed some characteris 
tics of Z sources near the peak of its outburst, which 
reached ^L-^dd (other bright NS-LMXB transients typi- 
cally only rea ch ^O.SLEdd). and atoll propertie s at lower 
luminosities (Wijnands & van der Klis 1999). While 



these two examples suggested that a single source can 
show atoll/Z source characteristics depending on lumi- 
nosity, and lent some credence to the belief that the NS- 
LMXB subclasses are entirely the result of differences in 
M, the actual evolution from Z tracks to atoll tracks had 
never been observed. 

Significant progress on the issue of NS-LMXB sub- 
classes was gained following the outburst of XTE J1701- 
462, a transient NS-LMXB that was discovered in Jan- 



uary 2006 (Rcmillard et al. 2006) and remained in out- 
burst for ~19 months. Initial observations with the 
Rossi X-ray Timing Explorer (RXTE) revealed that 
XTE J1701-462 was the first transient NS-LMXB to 
show all of the Z source characteristics in X-rays (H07). 
The radio properties of the so urce were consisten t with 
those of the Z sources as well ( Fender et al.| [2007). Ob- 
servations of type I X-ray bursts with radius expansion, 
near the end of the outburst, put the source at a dis- 



XTE J1701-462 



tance of 8. 8± 1.3 kpc and indicate that during the early 
phase of the outburst the sou rce had super-Eddington 
luminosities (Lin et al. 2009a). After the early phase 
of the outburst the source made a transition from Cyg- 
like to Sco-like Z source behavior as the source intensity 
dropped by a factor of 1.5 (H07), and observations near 
the end of the outburst, in July 2007, showed the source 
evolving f rom Sco-like Z source b ehavior to atoll source 
behavior (Homan et al. 2007a|c ). A full spe ctral anal- 
ysis of the entire outburst was performed by |Lin et al.| 
(2009b hereafter LRH09), who confirm that at the low- 
est luniinosities the spectral behavior of XTE J1701-462 
is consistent with that observed for the two atoll source 



transients Aql X-1 and 4U 1608-52 (Lin et al. 2007) 



XTE J 170 1-462 is the first source in which the evolution 
from Z to atoll source has been observed in detail, and 
its behavior strongly suggests that different levels of M 
are indeed sufficient to explain the different NS-LMXB 
subclasses. 

The outburst of XTE J1701-462 also allows for a study 
of the evolution of the X-ray properties along CD/HID 
tracks within the context of the different NS-LMXB sub- 
classes. This may shed new light on the nature of the 
Z and atoll branches, whether some or all are unique to 
each class, and, most importantly, what drives the mo- 
tion along them. As far as atoll sources are concerned, 
past observations of transient NS-LMXBs have suggested 
that the evolution along atoll source CD / HID tracks is 
mainly driven by var iations in M (see, e.g., van der Klis 
Lin et al.|2007 ) , increasing from the extreme island 



1994 



state to the upper banana branch. However, the afore- 



mentioned hyst eresis in spectral state transitions (Mac 
carone fc Coppi 2003) and parallel tracks observed m 



carone &: L^opp i zUUci) ana parallel tracks observed m 
HlDs (see, e.g., |l\fengez et al.|2001[|Di Salvo et al.|2003 l 



hint at a more complex relation between M and spectral 
state in atoll sources (Linares 2009), not unlike what 
is suggested by the spect ral evolution seen in transient 
black- hole X-ray binaries ( Homan fc Belloni|[2005j ) . 

It should be noted that there is ambiguity m some of 
these discussions about the definition of M. Some au- 
thors implicity define M as the mass transfer rate from 
the companion to the compact object Roche lobe, while 
others take the mass flow rate through the inner disk or 
the accretion rate onto the neutron star surface as the 
more relevant definition. Others explicitly distinguish 
various M components. In the presence of the various 
proposed fiows in the accretion/ejection process, such as 
disk and equatorial boundary layer fiows, spherical in- 
flows, magnetically dominated polar flows, possible disk 
winds and jets, clearly this is a complex issue to which 
we shall return in fJ4] 

Interpreting the evolution along the Z source tracks 
has been more difficult, owing to the fact that motion 
along the Z source tracks corresponds to relatively small 
(less than a factor of ~2) and non-monotonic intensity 

variations. For Z sources it was suggested that M in- 
creases from the horizontal branch to the fiaring branch 
(Hasinger et al. 1990} , although other scenarios have 
been put forwar d as well: e.g., M changing in the op- 
posite directi£nJChurc^et^^jj2006), or M not changing 



at all ( Homan et al.|20U2[ ). Based on their spectral anal 
ysis of the Z tracks, LKHi 



the Z tracks may occur at a nearly constant Af , with the 
branches being the result of different instabilities in the 
near-LEdd accretion fiow. 

More detailed analyses of the data presented in this 
work, dealing with different aspects of XTE J1701-462, 
can be found in a number of other pa pers: spectr al anal- 
ysis (LRH0 9), type-I X-ray bu rsts ( |Lin etal.| [2009a), 
kHz QPOs (Sanna et al. 2010), broad- band variability 
(Aresu et a L 2010, in prep.), and t he rapid decay into 
quiescence (Fridriksson et al. 2010). Here, we present 
an overview of the outburst of X'l'E J1701-462 that fo- 
cuses on aspects of the evolution of XTE J1701-462 in 
which the change from Z source to atoll source is most 
clearly seen: the CD/HID tracks, high-frequency QPOs, 
and broad-band variability. In particular, we find that 
selections of data groups based on low-energy count rate 
allow for a more detailed study of the evolution of the 
CD/HID tracks than in H07 and LRH09. This also per- 
mits a more precise comparison of the observed tracks 
for XTE J1701-462 with those found in the various NS- 
LMXB sub-classes. In Sj2]wc summarize our data set and 
analysis techniques. Our results are presented in f|3]aiid 
in SJ4] these results are interpreted and discussed within 
the framework of various scenarios for the role of M in 
driving the evolution between NS-LMXB subclasses and 
along the CD/HID tracks of NS-LMXBs. 

2. OBSERVATIONS AND DATA ANALYSIS 

We analyzed all 865 pointed RXTE pradt et al |1993| ) 
observations of XTE J1701-462 made between 200o Jan- 
uary 19 (MJD 53754) and 2007 Aug 29 (MJD 54341). 
Five of them were discarded from further analysis be- 
cause proportional counter unit 2 (PCU2) was not work- 
ing or because the length of usable data intervals was 
too short (<256s) for our variability analysis. The re- 
maining 860 observations had a total combined exposure 
time of ~2.73 Ms. For our analysis we only made use 
of data frora the Proportional Counter Array (PCA; Ja- 



hoda et al 

as described in 



i6|). 
■HO 



20061). We followed the same analysis steps 
07, the only difference being that in ad 



109 suggest that motion along 



dition to background corrections, the standard-2 data 
used for the light curves and CDs/HIDs in j ]3.1| were also 
corrected for dead time. All count rates and colors given 
in the text and used in the figures are for PCU2 only; 
for our power spectral analysis we made use of all ac- 
tive PCUs. Dates will be referred to as days since 2006 
January 19 (MJD 53754). 

3. RESULTS 

3.1. Light curves and color-color diagrams 

In Figure [2] we plot RXTE/PCA light curves of XTE 
J1701-462 in two energy bands: -2-2.9 keV and ~9.4- 
18.1 keV, corresponding to staiidard-2 channels 1-3 
and 20-40, respectively. The two light curves have a 
strikingly different appearance. Except for the first 30 
days of the outburst and occasional short drops in in- 
tensity between days 30 and 130 (see inset in Fig. [2Ji), 
the low-energy light curve shows little short-term vari- 
ability. Smooth long-term modulations with periods of 
—20-50 days are present, until the source starts its de- 
scent into quiescence. The high-energy light curve, on 
the other hand, shows strong short-term fiaring, which, 
as we discuss later, corresponds mostly to motion along 
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Figure 2. RXTE/PCA light curves of XTE J1701-462 in two energy bands; -2-2.9 keV (top) and ~9. 4-18.1 keV (bottom). The count 
rates are from the PCU2 detector and the time resolution of the light curves is 256 s. Error bars are smaller than the plot symbols. The 
inset in the top panel zooms in on the first few long-term cycles and the vertical lines near day 550 represent the start ti mes o f the three 
type I X-ray bursts that were detected. The count rate levels used for our CD/HID selection process, which is described in ^3.1[ are shown 
on the right-hand side of the top panel. Red data points are from observations (on the Z source horizontal branch) that were moved to a 
different CD/HID in this selection process. 



the Z source flaring branch. This flaring starts to weaken 
around the time at which the long-term low-energy mod- 
ulations end, and it abruptly subsides when the decay at 
low energies accelerates (day ~550). After day ^560, 
the count rates reach a very low, but non-zero, level of 
^2 ctss~^ per PCU (full energy band). This residual 
emission can be attribut ed to diffuse Galactic emission 
( Fridriksson et al. 2010), but has not been subtracted 



from our data, i'he times of the three type I X-r ay bursts 



that were detected near the end of the outburst ( Lin et al. 



2009a) are indicated by the short vertical lines m i" 'igure 

As shown by H07 and LRH09, XTE J1701-462 displays 
strong secular changes in its CD and HID. Creating a sin- 
gle CD/HID for the entire outburst results in a compli- 
cated set of overlapping tracks with varying shapes (see, 
e.g.. Figure 6 in LRH09). Time-based selections, such as 
the ones made by H07 and LRH09, also encounter strong 
secular motion associated with the 20-50 day modula- 
tions, which causes the shape and position of the tracks 
to change on a time scale of a few days. Fortunately, we 
have found that the low-energy count rate (2-2.9 keV; 
see Fig. [2^) appears to be a fairly good tracer of the sec- 
ular changes. For most of the outburst, the count rate 
in this band changes very little as the source traces out 
a particular Z track (although it does change along the 
atoll track we observed at the end of the outburst). Fur- 
thermore, the low-energy count rate can be used to make 
separate CDs and HIDs based on count rate intervals se- 
lected over most of the outburst. 

Eleven low-energy count rate intervals were defined to 



group the CD/HID data points; the maximum values of 
these intervals are shown as horizontal lines on the right- 
hand side of Figure [2^. Although the goal of these se- 
lections is to follow the secular evolution of the CD/HID 
tracks, there are two cases where changes in the low- 
energy count rate were mainly the result of motion along 
the CD/HID track. The first case concerns the period 
before day 29, during the Cyg-like Z source stage of the 
outburst, when the source shows large intensity swings in 
the low energy band, associated with motion along the 
Z. For this period we manually selected data to create 
two sets of CDs/HIDs (bringing the number of CD/HID 
selections to a total of thirteen) ; the first period consists 
of data from days 14.0-23.0, corresponding to the highest 
observed count rates, and the second combines data from 
days 0-14.0 and 23.0-29.0. The second case where we ex- 
perienced difficulties with the low-energy countrate selec- 
tions consists of several occasions when the source enters 
the horizontal branch after day 29. These instances are 
marked in red in Figure [2] and many are clearly visi- 
ble as drops in count rate near the maxima of the long- 
term modulations, as the result of which they end up in 
a lower count rate selection. To correct for this, these 
data points were grouped with their neighbouring non- 
horizontal branch observations (typically one count rate 
level higher). 

Twelve of the thirteen CDs and HIDs are shown in 
Figures [3] and |4j respectively, in order of decreasing low- 
energy count rate. The average total count rate of the 
tracks decreases (monotonically) in the same direction as 
well. The twelve CDs and HIDs will be referred to as se- 
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Figure 3. The evolution of the CD tracks of XTE J1701-462 as a function of low-energy count rate. The selection criteria for the tracks are 
given at the bottom of each panel, together with the total amount of exposure time. Diagrams C-L are based on selections of the low-energy 
(~2.0-2.9 keV) count rate R\ diagrams A and B arc based on manual selections (before day 29). Soft color and hard color are defined 
as the ratio of count rates (extracted from standard-2 data) in the ~4. 0—7.3 keV (channels 7-14) and '~2. 4-4.0 keV (ch. 3-6) bands, and 
the ~9.8— 18.1 keV (ch. 21-40) and ~7.3-9.8 kcV (ch. 15-20) bands, respectively. Data points represent 64 s averages for high count rates, 
increasing to 1024 s averages for the lowest count rates. The blue points indicate observations in which kHz QPOs were detected; the red 
points represent data from observations on th e Z source horizontal branch that were moved to the selections corresponding to neighboring 
non-horizontal branch observations (see gSTTJ. The dashed gray line indicates the approximate path followed by the normal/flaring branch 
vertex during the transition from Cyg-like L to Sco-like Z and, finally, atoll source behavior. 



lections A-L throughout the remainder of the paper. Ex- 
cept for observations with low-energy count rates below 
0.4 ctss~^ (i.e., the thirteenth count rate interval), which 
have data with large errors bars, these selections repre- 
sent our entire RXTE data set of XTE J1701-462. The 
count rate intervals used for each CD/HID are shown 
near the bottom of each panel, except for selections A 
and B, which were based on manual selections and had 
peak low-energy count rates of ^^100-165 ctss~^. The 
total amount of time in each selection is shown as well. 
The horizontal branch observations marked with red data 
points in the light curves in Figure [2] are shown in red 
in Figure l3] as well. Compared to the CDs shown in 
LRH09, which were based on broad time selections, the 
set of CDs shown in Figure |3] reveals a more detailed evo- 
lution of the CD tracks. We note that the ordering of the 
panels in Figures |3] and |4] does not reflect the evolution 
of the source as a function of time; as can be seen in Fig- 
ure [2)!i, XTE J1701-462 moved back and forth between 
different selections (C-H) for most of the outburst. 

While our low-energy count rate intervals result in rel- 
atively narrow CD tracks, there is still considerable sec- 
ular motion within some of the HID tracks, making it 
harder to follow the evolution of these tracks than in the 
CDs. For this reason we created subset HIDs from nar- 
rower count rate selections within selections D-K, which 
are shown in red in Figure |4J 

Significant evolution can oe seen in the shape of the 
CD/HID tracks. It has already been reported that the 
CD/HID tracks of XTE J1701-462 evolved from Cyg-like 
Z tracks, to Sco-like Z tracks, an d finally into atoll-like 
tracks (II07, [Homan et "al}j2007c^ LRII09), which can 
clearly be seen in !< igures |5[ and |4[ Typical Z and atoll 
branches have been marked in a few of the panels of 
Figure [3] In the following we take a more in-depth look 
at the evolution of the CD/HID tracks and also compare 
the CD/HID tracks of XTE J1701-462 to those found in 
other NS-LMXBs. 

3.1.1. Comparison with CDs/HIDs of other NS-LMXBs 



To put the evolution of the CD/HID tracks of XTE 
J1701-462 into a broader context, we compare the var- 
ious CD/HID tracks observed in XTE J1701-462 with 
those seen in other NS-LMXBs. For this comparison we 
mainly refer to work that is based on RXTE/PGA data, 
as comparison with data from other satellites is com- 
plicated by differences in energy bands and instrument 
responses. We also refer to the CDs/HIDs presented in 
Figure 111 which were taken from Fridriksson et al. (2010, 
in prepT). 

• Selections A & B: The overall shapes of the 
CD/HID tracks in these two selections resemble 
those see n for Cyg X-2 in its high l uminosity state 
(Wijnan ds fc van de r Klis 2001 ; Muno et al.||2002' 



Piraino et al. | |2002| |Leret~ai:[|20U8| ). They~aIso 
show similarities to the tra cks observe d in the Cyg- 
like Z sources GX 340-^0 ( [Jonker et al. 2000), GX 
5-1 (Jonker et al. 2002 see also I'ig. Ilk), how- 
ever, with some notable differences in the prop- 
erties of the flaring branch. The flaring branch 
in the CD of these two sources is usually oriented 
horizontally, rather than pointing downward like in 
XTE J1701-462, although downward-pointing ex- 
tensions to the flaring branch have been observed 
for these two sources (the so-called 'dipping' flar- 
ing branch). Moreover, in GX 5-1 and GX 340+0 
the count rate increases when they move onto the 
flaring branch (at least up to the 'dipping' flaring 
branch), whereas the count rate decreases strongly 
on the flaring branch in selections A and B. Such 
dramatic drops in intensity on the flaring branch 
have only been seen in Cyg X-2 during its high lu- 
minosity state. We note that the track of selection 
A is incomplete, as evidenced by the visible gaps. 

• Selection C: This selection contains relatively little 
data; the track is likely to be incomplete and there- 
fore difficult to compare with other sources. How- 
ever, the properties of the flaring branch in the CD 
(short length and horizontal orientation) and HID 
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Figure 4. The evolution of the HID tracks of XTE J1701-462 as a function of low-energy count rate. Black data correspond to the 
same data as shown in Figure Is] whereas the red data represent narrower count rate intervals with reduced effects of secular motion. The 
selection criteria for the tracks are given at the bottom of each panel, together with the total amount of exposure time (only for the entire 
selection). Diagrams C-L are based on selections of the low-energy (~2.0-2.9 keV) count rate R; diagrams A and B are based on manual 
selections (before day 29). For the intensity we used the total PCU2 count rate; for the definition of hard color see the caption of Figure p] 
Data points represent 64 s averages for high count rates, increasing to 1024 s averages for the lowest count rates. The blue points indicate 
observations in which kHz QPOs were detected. The dashed gray line indicates the approximate path followed by the normal/flaring branch 
vertex during the transition from Cyg-likc Z to Sco-like Z and, finally, atoll source behavior. 



(count rate increase) suggest that it might be part 
of a track more similar to those of GX 5-1 and GX 
340+0 than the ones seen in selections A and B. 

• Selections D to H: These selections show CD tracks 
that are similar to those of the Sco-like Z sources 



GX 17-F 2 dHOTmn et al.|[2002l see also Fig. ^) 
ScoX-1 ( van der Khs et aT 1199 6' 'Bradshaw cT~al. 
2003D, and"GX 849-^2 ([O'JNeill et al. 2002; Agrawal 



&: Sreekumar 2003 see also b'ig. fl|:) ; they also show 



GX l'6+l ^Muno et al.) |2 
The latter is otten classini 



simila rities to those of Cyg X-2 in its low luminosity 
state dKuu lkers et al.||l 996; Muno et al. 200 2| and 

iSchncrr ct al.||2003[ ). 
ed as a bright atol l source, 
but also dis plays Z source characteristics ( [Homan 
et al. [2004 1 ■ Note that in the HIDs of selections 
D and E the horizontal branch is still relatively 
long and close to being horizontal, similar to the 
Cyg-like Z sources, but unlike what is observed for 
the Sco-like Z sources. The properties of the hor- 
izontal branch become more Sco-like in the HID 
of selection F, where it turns upward and short- 
ens. Large intensity increases associated with the 
flaring branch, which are typical for the Sco-like Z 
sources, only start to become visible in selections 
E and F. 

• Selection I: In this selection only hints of the nor- 
mal branch can be seen in the CD and HID, in addi- 
tion to a full-fledged flaring branch. T he EXOSAT 
observations of GX 349-1-2 analyzed by |Hasinger fc| 
van der Klis ( 1989 1 reveal a similarly shaped track. 
Overall, the track appears to be intermediate to the 
Z source flaring branch and the banana branches 
displayed by bright atoll sources such as GX 9-1-1 
and GX 9+9 (see below). 

• Selections J & K: No normal branch is seen in 
these selections and the curvature in the remain- 
ing branch in the CDs, while still retaining some 



characteristics of the flaring branch (particularly 
in selection J), is also reminiscent of the lower and 
upper banana branches displayed by brigh t atoll 
sources such as GX 9+1 and GX 9+9 (Hasinger & 
van der Klis|[T989 see also Fig. IJi) and to a lesser 
extent GX 3+1. Curvature of the track is stronger 
in selection K, where the lower banana branch be- 
comes more prominent. The intensity shifts shown 
by the (upper) banana branch in the HIDs, which 
are the result of secular motion, are also seen for 
the bright atoll sources; in particular, the HID of 
selection K is similar to that of the atoll source 4U 
1735-44, which also shows secular motion of the 
upper banana branch (Fig.[T]3). 

• Selection L: This selection shows a combination of 
branches (lower-banana, island state, and extreme 
island state) that is frequently seen in persistent 
atoll sources such as 4U 1636-53 and 4U 1705-44 

see also 



Homan et al. 2009 



(jBelloni et aL||2007| 

Fig. [ijjg) and transie nt atoll sources such as Aql^ 



X-1 and 4U 1608-52 (Lin et al. 2007 Gladstone 



et al. 2007). For a better comparison with these 



sources we show the comb ined CDs and HIDs of 
selections J-L in Figure [5J [Gladstone et al. ( 2007 1 
found considerable variation in the orientation of 
the island state among the atoll sources they stud- 
ied. In XTE J1701-462 the island state is vertical 
in the CD and pointing to the upper left in the HID, 
similar to 4U 1735-44 (Fig. ^) 



( [Gladstone et al.[[2007[ ) 



and 4U 0614+09 



3.1.2. Evolution of tracks and branches 

The above comparison with other NS-LMXBs con- 
firms our earlier reports that the CD/HID tracks of XTE 
J1701-462 evolved from Cyg-like Z to Sco-like Z, and fi- 
nally became atoll-like. What follows is a brief overview 
of the most significant aspects in the evolution of the 
CD/HID tracks and branches. 
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Figure 5. The combined CD (a), HID (b), and 9.4-18.1 keV light 
curve (c) of selections J— L The red-colored data points show tliat 
the part of the atoll source track that has historically been classified 
as the "upper banana" corresponds to the last episodes of flaring 
in the Z source phase. This flaring abruptly subsides below total 
count rates of ~300 ctss~^ (b), corresponding to a high energy 
count rate of ~40 ctss"'^ (c). The numbered circles in the CD 
and HID show the location of the three type I X-ray bursts that 
were detected. The times of the three bursts are indicated by the 
three vertical lines in the light curve; the bursts themselves are not 
shown. 



The horizontal branch in the high-intensity and Cyg- 
hke selections A and B shows large intensity swings that 
only results in minor color changes. Selections A and B 
are also the only two selections to show an extension to 
the horizontal branch that previously has been referred 
to as the horizontal branch 'upturn' (H07, LRH09). In 
the CD this upturn appears to be more of a separate 
branch than in the HID, similar to what is seen in GX 
5-1 an d in the EXOS AT data of Cyg X-2 ( |Kuulkers et al 
19961. The horizontal branch is not observed in selection 
C, probably due to the short total exposure time. In the 
CDs of selections D-F the orientation of the horizontal 
branch is rotated in a clockwise direction with respect 
to selections A and B, while the intensity swings have 
become less pronounced. The upturn is no longer present 
in these selections, and the horizontal branch itself is no 
longer visible in selection G. 

The normal branch gradually becomes shorter in the 
CDs and HIDs from selection A to H, while it rotates 
slightly in a counter-clockwise direction. As in the case 
of the horizontal branch, the shortening of the nor- 
mal branch in the HID means that the intensity swings 
become less pronounced. The normal branch is only 
marginally detected in selection I. As was already pointed 
out by H07, the vertex between the normal and flaring 
branches moved along a more or less straight line in the 
CD during the first part of the outburst, as the result of 
the secular changes. As can be seen in Figure [3J where 
we have represented the path of the vertex with a dashed 
line, this also appears to be true for the remainder of the 
Z source phase of the outburst. LRH09 showed that this 
vertex also followed a straight line in their HID. With 
our color definition (i.e., chosen energy bands) the path 



followed by the vertex in our HID shows some deviations 
from a straight line (especially in the three highest count 
rate selections), but it is still rather close. The dashed 
line in Figure |4] represents the approximate path of the 
vertex; it is a straight line, but it appears curved because 
of the logarithmic scale on the horizontal axis. 

The flaring branch undergoes large changes in both the 
CD and HID and is crucial in linking Z and atoll source 
behavior. In selections A and B motion onto the flar- 
ing branch corresponds to drops in the intensity, and in 
the HID the flaring branch appears to be an extension 
of the normal branch. In the CD of selection A the flar- 
ing branch points downward, but it rotates in a counter- 
clockwise direction in the lower count rate selections. 
Starting in selection C, motion onto the flaring branch 
corresponds to count rate increases, which become more 
pronounced in selections E and F. This flaring eventually 
starts to become weaker (see also the high-energy light 
curve in Figure Pp), but it remains present through se- 
lection K. While selections J and K can be classified as 



atoll-like (^3.1.1), it is important to note that the flar- 
ing that has long been regarded as being characteristic 
of the Sco-like Z sources is still present in these lower 
count rate intervals. In the framework of atoll sources, 
the branches traced out as a result of this flaring would 
be referred to as the atoll upper-banana branch, but as 
Figures [2J [3 and 4] show, there is no sharp transition be- 
tween Z and atollbehavior. To better describe the rapid 
evolution during the last phase of the outburst and em- 
phasize the relation between the Z source flaring and the 
atoll upper-banana branch, we show the combined CD, 
HID, and light curve of selections J, K, and L in Figure 
[5] resulting in a 'complete' atoll track. Points on the 
upper-banana branch of the atoll in the CD were col- 
ored red; from the HID and light curve it can be seen 
that these upper-banana branch data points correspond 
to the last few flares that were observed in XTE J1701- 
462 (compare also with Figs.llJD and 111). This is the first 
time that the Z source flaringbrancnnas been observed 
to evolve into the atoll source up per banana branch and 
it c onfirms earlier suspicions (Hasinger & van der Klis 
|1989j ) that the two branches might be one and the same 
phenomenon. 

The lower banana branch starts to form around the 
time of the last flares in selections J and K. As the total 
count rate drops by a factor of ^3 from selection K to se- 
lection L, the source moves away from the diagonal line 
followed by (the secular motion of) the normal/flaring 
branch vertex in selections A to I (dashed lines in Fig. 
[3]). However, in the HID the lower banana branch con- 
tinues along the path traced out by the secular motion of 
this vertex (dashed lines in Fig.l4|, as already pointed out 
by LRH09. The transition to the atoll island state oc- 
curs in selection L during a rapid decrease in count rate, 
and is marked by a sudden increase in hard color. It is 
followed by the extreme-island state, in which the hard 
color becomes more or less constant (although poorly 
constrained because of the low count rates) . 

3.2. Summary of CD/HID evolution 

Figures [3] and H] show that XTE J1701-462 is the flrst 
NS-LMXB in which we can study the entire evolution 
from Z source to atoll source behavior in a single source. 
Cyg-like Z source behavior, characterized by large inten- 
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Figure 6. Power spectra from the lower-banana (a) and upper- 
banana branch (right) of the GX-atoU-like selections I and J, shown 
with their best fits. The feature around 160 Hz in panel (a) is 
some what suggestive of a "hectohertz QPO" J van Straaten et al.| 
|2002[ |, but it is only detected at a 1.8 a level. 

sity swings on the horizontal and normal branches, the 
presence of a dipping flaring branch, and an upturn from 
the horizontal branch, is observed at the highest count 
rates. The transition to Sco-like Z source behavior is 
rapid (~5 days) with maximum total count rates drop- 
ping by a factor of ^1.5; this transition mostly involves 
changes in the orientation of the branches in the CD 
and HID, as well as the disappearance of the horizontal 
branch upturn. The transition from Sco-like Z behav- 
ior to atoll source behavior is much slower (~months). 
Figure [5] shows that the last instances of the Z source 
flaring branch, at low overall count rates, can also be 
classified as the atoll source upper banana branch, pro- 
viding an important link between the two subclasses. 
The Sco-like Z to atoll transition is marked by the dis- 
appearance of the horizontal branch, normal branch, 
and flaring/upper-banana branch, as the maximum total 
count rates dropped by a factor of ~3. In terms of the 
selections A-L this is a systematic and monotonic pro- 
cess; of course, during its initial rise and in subsequent 
decay the source moves back and forth between some of 
these selections. The time between the last detection of 
the flaring branch/upper banana branch to the transi- 
tion from the lower-banana branch to the island state 
(and later the extreme island state) is ~7 days, with the 
count rate dropping by another factor of ^3. The evolu- 
tion during the atoll phase presents an important depar- 
ture from the evolution seen in the Z phase of the out- 
burst. While in the Sco-like Z source phase the change 
in low-energy count rate results in gradual changes in 
the CD/HID tracks (except for occasional visits to the 
horizontal branch), in the atoll phase (after the flaring 
subsides) the low-energy count rate changes result in mo- 
tion along the CD/HID track. In the HID, motion along 
the lower-banana branch even appears to be a continua- 
tion of the gradual motion of the Z source normal/flaring 
branch vertex (this is not the case in the CD, however). 

3.3. Low-frequency variability 

A low-frequency variability study (in the 6.9-60 keV 
band) of the Cyg-like tracks and some of the Sco-like 
tracks was already presented in H07. Two types of low- 
frequency QPOs were observed in those tracks, the hor- 
izontal branch oscillations (~10~60 Hz) and the normal 
branch oscillations (~7-9 Hz). These QPOs appeared on 
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Figure 7. The island state power spectrum of atoll-like selection 
L. The solid line shows the best fit (consisting of three Lorentzians) 
to the power spectrum. 



the appropriate branches, confirming our branch identi- 
fication and our interpretation of the highest count rate 
tracks as Z-like tracks. We did not perform a similarly 
detailed analysis for the remainder of the outburst, but 
a visual inspection of the power spectra indicates that 
both types of low-frequency QPOs disappear as the over- 
all count rate decreases. The normal branch oscillations 
disappear within selection E, where they are seen in the 
high count rate HID tracks, but not the lower count rate 
ones. Indications of horizontal branch oscillations are 
still seen in selections F and G, but not in the lower 
count rate selections H-L. 

In Figure [6] we show examples of 2-60 keV lower- 
banana (a) and upper-banana branch (b) power spec- 
tra from the combined bright-atoll-like selections I and 
J. The shapes are consistent with the power spectra 
of ot her atolls in these spectral states (van der Klis 
2006 1 . The lower-banana power spectrum shows a broad 
(Q=0.7±0.1) peak at 13.9±0.7 Hz, on top of a power- 



law shaped (ex v 



-1.41±0.04 



) noise continuum. The over- 



all variability is weak, with an integrated power rms 
(1/16-100 Hz) of ~3% (2-60 keV). The upper-banana 
power spectra could be fitted with a single power law 
(oc ;y-i-48±o.03) ^^^ tjjg 1/16-100 Hz rms was ~2.4%. 

A similar study for the atoll-like selections K and L 
suffers from low count rates, and high quality broadband 
power spectra are difficult to produce. The shape of 
the power spectrum of the lower-banana branch of se- 
lection L (from observations with total count rate above 
65 counts/s/PCU) is consistent with that of the lower- 
banana branch of selections I and J; we measure an rms 
amplitude of 5.0ib0.3% . which is also consistent with that 
seen in atoll sources (van der Klis 2006). The island 
state power spectrum was averaged from observations 
with hard color 0.9-1.1 and total count rate above 35 
counts/s/PCU. The resulting power spectrum is shown 
in Figure [T] It could be fit with a model consisting of 
three broad (0 < Q < 1.8) Lorentzians, with frequencies 
of 27±4, 155±19, and 720±50 Hz. The 720 Hz feature is 
probably a (broadened) lower or upper kHz QPO (see be- 
low), whereas the 155 Hz feature could be a hecto-hertz 
QPO, a component that i s frequently seen in this state 
( |van Straaten et al.|[2002l [2003| ). The broad 27 Hz fea- 
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ture likely corresponds to a low-frequency band-limited 
noise component. We measure a 1/16-100 Hz rms ampli- 
tude (2-60 keV) of 16.3±0.7%, which combined with the 
overall shape of the power spectrum is consistent with 
the island state of other atoll sources. For the extreme 
island state we averaged power spectra from observations 
with total count rates between 4 and 19 counts/s/PCU 
(with hard color ^1.1-1.3). We measure an rms ampli- 
tude of 29±2%, also similar to what is found in the ex- 
treme island state of other atoll sources. For reference, 
the highest rms for the Z-like selections was measured 
on the horizontal branch upturn in selections A and B: 
15.0±0.2%. 

3.4. Kilohertz QPOs 

Following H07, the 2-60 keV and 6.9-60 keV power 
spectra of all individual observations were visually in- 
spected for high-frequency (kHz) QPOs. In addition 
to the (marginal) detections reported in H07, we find 
two pairs of kHz QPOs and six single kHz QPOs. The 
observation IDs of all ten > 3a kHz QPO detections 
(pairs and singles) are listed in Table [IT together with 
their CD/HID selection and QPO frequencies; frequen- 
cies are from fits to the power spectra of the 6.9-60 keV 
energy band. The locations of these observations in the 
CDs/HIDs are shown in blue in Figures Is] and [4] One 
new pair was detected on the horizontal Drancn of se- 
lection D and the other was detected in two consecutive 
observations (made about a day apart) on the horizon- 
tal branch of the Sco-like track of selection F; the power 
spectrum of the second pair is shown in Figure [Sk. The 
QPO properties are consistent with those of trie pairs 
reported in H07, which were detected on the horizontal 
branch or close to the horizontal/normal branch vertex 
of the Sco-like CD/HID tracks of selections D, E, and 
F. The kHz QPO pairs have peak separations between 
274±11 Hz and 319±7 Hz, and the average peak separa- 
tion is ~285 Hz, consistent with wh at is found in other 
NS-LMXBs ( [Mendez fc Belloni|2007"l ). The six new single 
kHz QPO detections were all rnade on the lower-banana 
of the atoll-like CD/HID track of selection L. The 720 
Hz fe ature in the combined island state power spectrum 
(j |3.3[ ) was not detected significantly in individual obser- 
vations. Besides the fact that these single peaks were de- 
tected in an entirely different part of the CD, their prop- 
erties also differ considerably from the peaks detected 
earlier in the Sco-like CD tracks. As can be seen from 
the representative example in FigurelSb, these peaks have 
rms amplitudes that are two to three times higher than 
the peaks in the Z like tracks and also have much higher 
Q-values, consiste nt with the (lower) kHz QPO s found 
in atoll sources (M endez|2006 van der Klis|2006|. Their 
frequencies ranged from ~650 Hz to ~'850 Hz, and in 
dynamical power spectra of two individual observations 
(with lengths of ^^8-15 ks) the frequ e ncy w as found to 
drift by up to ~80 Hz. Sanna et al. (2010), combining 
all kHz QPO data froni the lower banana branch and 
employin g a shift-and- add technique to align the drifting 



peaks ( jMendez et al. 1998), report the marginal detec 

tion of a second peak, with a peak separation on 258±13 

Hz. Finally, we note that averaging power spectra based 
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Figure 8. Two examples of kHz QPO detections in the 6.9-60 keV 
band. On the left we show a pair of kHz QPOs from the horizontal 
branch of selection F (observations 92405-01-40-04 and 92405-01- 
40-05 combined). On the right we show a single kHz QPO from the 
lower banana branch of selection L (observation 93703-01-02-05). 
Fit parameters for all QPOs are shown. Notice the difference in 
vertical scale between the two panels. 



Table 1 

Kilohertz QPO detections in XTE J1701-462 



Observation ID 


Selection 


u (Hz) 


91442-01-07-09 


D 


641±5, 915±10 


92405-01-01-02 


F 


613±4, 932±6 


92405-01-01-04 


E 


761±6 


92405-01-02-03 


D 


622±9, 911±12 


92405-01-03-05 


F 


617±11, 908±9 


92405-01-40-04° 


F 


644±5, 923±6 


93703-01-02-04 


L 


852±2 


93703-01-02-11'' 


L 


~805 


93703-01-02-05 


L 


834±1 


93703-01-02-08 


L 


837±2 


93703-01-03-00'= 


L 


~720 


93703-01-03-02 


L 


650±8 



on selections in the CD/HID (e.g., Homan et al. 20021 



revealed additional indications for kHz QPOs on the hor 
izontal and normal branches of the Sco-like selections D- 



" Fitted together with power spectrum of observation 92405-01-40- 

05 

*" QPO moved by ~60 Hz during observation 

*= QPO moved by ~80 Hz during observation 

F; however, a full study of this kind is beyond the scope 
of this paper. 

4. DISCUSSION 

4.1. Identification of Z and atoll phases 

In ^ we presented an overview of the 2006-2007 out- 
burst of the transient NS-LMXB XTE J1701-462. By 
comparing the CD/HID tracks of XTE J1701-462 with 
other sources, we showed that the source made a rapid 
transition from Cyg-like to Sco-like Z source tracks, fol- 
lowed by a slower transition to various types of atoll 
source tracks. The transitions between Z and atoll and 
further subclass divisions are not monotonic but closely 
track the 2-2.9 keV count rate; we note that the spec- 
tral fits by LRH09 suggest that for most branches the 
count rate in this energy band is dominated by the con- 
tribution from the accretion disk component. During the 
outburst, the count 2-2.9 keV rate initially goes down 
rapidly, which is linked to the rapid Cyg-like to Sco- 
like transition; later it goes up and down superimposed 
on the gradual decay, and this is associated with the 
source moving back and forth between subclasses. We 
interpret the 2-2.9 keV count rate as a proxy for the 
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mass accretion rate component M, possibly that from 
outer disk to the inner parts of the accretion flow, that 
determines subclass. XTE J1701-462 is the first NS- 
LMXB in which such transitions between subclasses in 
the CD/HID have been followed in detail (see also H07, 
LRH09). The properties of the low- and hig h-fr eque ncy 



source upper banana branch. This flaring rapidly be- 
comes weaker in the atoll phase, where it ev entual ly sub- 
sides around a luminosity of ~0.15LEdd (?3.1.2 lumi- 



variability observed in XTE J1701-462 (|3^and|3^ see 
also H07 and LRH09) broadly support our identification 
of the various Z-like and atoll-like tracks, although there 
are a few minor differences; w e exp ected to see kHz QPOs 
in the Cyg-like Z tracks (see [4.4) and normal branch os- 
cillations on the normal branch of selection F, but none 
were detected. 

Three type I X-ray bursts were detected in the e ntire 
RXTE dataset of XTE J1701-462 ( |Lin et al.|2009a| ). We 
have marked the times of the three bursts in the light 
curves in Figure [2] and also marked their location in the 
CD and HID in Figure [5] As can be seen, the bursts oc- 
curred late in the outburst, with two bursts occurring on 
the lower banana branch and one occurring on the atoll 
upper banana or Z source flaring branch (see { 3.1.2 1. The 
fact that the only three detections of type I X-ray bursts 
were made in the atoll phase of the outburst is in line 
with the fact th at atoll sources are mo re prolific bursters 
than Z sources (Galloway et al. 2008), and further sup- 
ports our ide ntification of the fi nal phase of the outburst 
as atoll-like. Lin et al. ( 2009a[ ) also derive a distance es- 
timate of 8.8±1.3 kpc from two radius expansion bursts 
(the latter two of the three), implying that the persistent 
luminosity just before the bursts was ^10% iEdd: also 
consistent with atoll sources. 

A spectral analysis of the complete RXTE data set of 
XTE J1701~462 is presented in LRH09. They find that 
near the end of the outburst, below luminosities of ^15% 
^Edd, the spectral behavior of XTE J1701-462 was sim- 
ilar to th at of the atoll-ty pe transients Aql X-1 and 4U 
1608-52 ( |Lin et al.|[2007J ). Notably, both the disk and 
boundary layer components closely follow L oc T on the 
atoll lower banana branch, implying a relatively constant 
emission area for both components as M changes. The 
spectral results of the Z source stage could not be com- 
pared to other Z sources, as the spectra of the latter have 
not yet been fitted with the model used by LRH09. 

4.2. The nature of NS-LMXB subclasses 

For each of the CD/HID tracks observed in XTE 
J1701-462 (Figs. [3) |4[ and [S]) similarly shaped track s 
can be found among the persistent NS-LMXBs (S 3.1.1). 



Moreover, except for sources that suffer from high inter- 
nal absorption (i.e., dipping and/or eclipsing sources), we 
are not aware of sources that show CD/HID tracks that 
are wildly different from the ones seen in XTE J1701- 
462. This seems to suggest that the evolution of the 
tracks seen in XTE J1701-462 is representative of the 
entire class of NS-LMXBs, which we will further explore 
in an upcoming paper (Fridriksson et al. 2010, in prep.). 
The gradual transition from Sco-like Z to atoll behavior 
observed in XTE J1701-462 argues against a sharply de- 
fined boundary between these two NS-LMXB subclasses. 
In particular, the flaring that is apparent in the high- 
energy light curves of XTE J1701-462 persists across 
the Z/atoll periods of the outburst, as it is responsi- 
ble for both the Z source flaring branch and the atoll 



nosity taken from LRH09). In addition to the disappear 
ance of the flaring, there are two additional points in the 
evolution of XTE J1701~462 at which significant changes 
occur; one is the transition from Cyg-like to Sco-like Z be- 
havior around LEdd, and the other is the switch between 
thermal- and non-thermal-dominated spectra (lower ba- 
nana/island state transition) around 0.02 Z^Edd (LRH09). 
While these points could serve as 'natural' boundaries be- 
tween different types of NS-LMXB subclasses, there are 
a number of persistent NS-LMXBs that regularly cross 
one or more of these boundaries, so their usefulness for 
(re-)defining NS-LMXB subclasses is therefore somewhat 
limited. 

As already mentioned in fJTl the disparity in the proper- 
ties of the Z and atoll subclasses has been a long-standing 
problem: are differences in mass accretion rate alone 
sufficient to explain the differences in the phenomenol- 
ogy of the various subclasses, or are additional param- 
eters needed? Differences in the neutron-star magnetic 
field have been suggested to be (partially) res ponsible for 
the differences bet ween Z and atoll sources (Hasinger & 
van der Klis||1989|) and betw een Cyg-like and Sco-like Z 
sources jPsaitis e t al.||1995|). Diamagnet ic screening by 



accreted matter ( Cumming et al.||2001 [) appears to be 
the only feasible way to affect (i.e. reduce) the external 
magnetic field of a neutron star on a short time scale. It 
cannot explain the existence of two classes of Z sources in 
terms of different magnetic field strengths, since fields in 
these sources should be completely buried as the result of 
their near-Eddington accretion rates. E xplaining the Z 
to atoll transition is problematic as well. [Hasinger fc van| 
der Klis^ (1989) suggested that atoll sources have a higher 
magnetic field than Z sources and to explain the observed 
Z to atoll transition, the underlying magnetic field needs 
to emerge again. However, the time scale on which this 
is expected to happen is much longer (100-1000 yr) than 
the observed time scale for the transition from Z to atoll 
in XTE J1701-462. We therefore conclude that magnetic 
field strength is probably not responsible for the NSXBs 
subclasses. 

Differences in viewing angle have been proposed to ex- 
plain the exi stence of two types of Z sources, Cyg-like 



and Sco-li ke (Kuulkers et al. 1994 
Klis 19951 and references therein 



Kuulkers & van der 
Changes in the bi- 
nary inclination are not expected on the time scale of 
weeks to months, so these can be ruled out as a source 
for the observed transitions between subclasses in XTE 
J1701-462. However, a tilted precessing disk, such as 
thoug ht to be present in, e.g., Cyg X-2 (Vrtilek et al. 



19881, can lead to changes in the viewing geometry on 
a tune scale of weeks to months. Combined with the 
effects of anisotropics and/or varying amounts of obscu- 
ration this could lead to changes in the morphology of 
CD/HID tracks and branches. However, such a model 
would require a complex evolution of the viewing angle 
to account for the observed sequence (in time) of the var- 
ious sub-classes. Moreover, it cannot account for the sys- 
tematic disappearance of the various Z source branches 
as is observed during the transition from Z to atoll. We 
conclude that viewing angle is probably also not respon- 
sible for the NSXB subclasses. We note that parameters 
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such as neutron star mass and spin frequency should not 
change substantially during an outburst, as the amount 
of matter that is accreted is negligible in comparison to 
the neutron star mass. 

Given the above, and taking also into consideration the 
transient nature of XTE J1701-462, we conclude that 
transitions between the NS-LMXB subclasses in XTE 
J1701-462 are probably entirely driven by changes in 
M . Although we cannot completely rule out that the 
persistent Z and atoll sources differ in, e.g., neutron star 
mass/spin/magnetic field, or viewing angle, our results 
indicate that these parameters are relatively unimpor- 
tant to the phenomenology of XTE J1701-462, and by 
extension, to that of the Z and atoll sources. Such pa- 
rameters may, however, explain some of the differences 
between tracks in the CD/HIDs of sources with similar 
M. 

The discussion of the role of the various possible in- 
flow/outflow components of the accretion process in de- 
termining NS-LMXB subclasses and motion along the 
CD/HID tracks is of course hampered by the lack of reli- 
able independent measures of the relevant inflow/outflow 
rates, so it is difficult to draw firm conclusions, even 
from the rich phenomenology observed in XTE J 1701- 
462. However, as we discussed above, for most of the 
outburst of XTE J1701-462 the low-energy count rate 
appears to be a fairly good tracer of the M that deter- 
mines the changes in the CD/HID tracks. It is plausible 
that this M scales monotonically with some overall rep- 
resentative average PCU2 (^-^2-60 keV) count rate level 
of the CD/HID tracks and their average bolometric lumi- 
nosities (LRH09). Under certain assumptions this would 
imply that Cyg-like Z source behavior corresponds to 
the highest overall energy release and perhaps M val- 
ues, followed by Sco-like Z source behavior and then atoll 
source behavior. The above also implies that the hori- 
zontal branch and its upturn are phenomena related to 
the highest M values, although, as we discuss below, this 
does not necessarily mean that along a single Z track 
these parts of the track represent the highest values of 
the accretion rate relevant to motion along the Z. 

4.3. Spectral evolution along Z and atoll tracks 

The above interpretation in which an M level deter- 
mines the type of CD/HID track of NS LMXBs has 
implications for proposed models for the role of mass 
flow rates in motion (i.e. spectral evolution) along the 
Z tracks. Such a picture is obviously problematic for 
models in which the same M level that determines the 
type of CD/HID also directly determines the position 
along the Z, with this M cither increasing monotoni- 
cally from tliejMri^ont^to_^ flaring branch (Hasingcr 
et al.||199"0 Vrtilek et al.||1990 1 , as is still often assumed, 
or with M changing in the opposi t e dire ction , as was 
recently pro posed by Church et al. (20061 and Jackson 
et al. (2009). However, we cannot rule out that the ac 



cretion rate through some part of the inner accretion flow 
changes monotonically along the Z. 

A suggested solution in which changes in a mass flow 
rate can explain both the changes in the type of the 
CD/HID tracks as well as motion along the Z track was 
discussed in H07. Their solution was based on a model by 



van der Klis (2001 ) for the so-call ed 'parallel tra.cks' p he- 
nomenon in INS-LMXBs (see, e.g., Mendez et al.|199"8t for 
further information on this), which assumes that there 
exists both a prompt and a flltered response to changes in 
the mass accretion rate Md through the inner edge of the 
accretion disk. For a detailed discussion of the flltered 
response scenario and how it relates to motion along the 
Z we refer to H07, but we point out an opportunity and 
a problem here. The model predicts that branches would 
be preferentially suppressed if there is a prolonged mono- 
tonic decay in M^, which is somewhat reminiscent of the 
systematic disappearance of the Z source horizontal and 
normal branches. However, the observed long term mod- 
ulation in 2-2.9 keV flux is not monotonic, and if it can 
be identifled with Md it should alternatingly suppress the 
horizontal and the flaring branch, which is not observed. 
An alternative scenario, in which motion alon g the Z 
track is not governed by M, was suggested by 



Homan 



et al. (2002). In such a scenario the parameter(s) re^ 
sponsible tor spectral and variability properties, e.g., the 
inner disk radius, varies (vary) independently from M. 
The spectral fits by LRH09 suggest that motion along 
the horizontal branch and the Sco-like flaring branch is 
consistent with occurring at a nearly constant M, al- 
though this is not as clear for the normal branch and 
Cyg-like flaring branch. In any case, the spectral flts sug- 
gest that M variations along individual Z tracks are very 
small. Based on this, LRH09 suggest that the different Z 
branches may be the result of various instabilities in the 
accretion flow, which could be related to an increase in 
the radiation pressure in different parts of the flow (inner 
disk, boundary layer). This scenario avoids the difficul- 
ties of Al having to explain both gradual changes in the 
shape of the Z tracks and motion along the Z track at the 
same time. The observations of XTE J1701-462 require 
that the occurrence of the assumed (and unknown) insta- 
bilities themselves depends on the M level (i.e., the insta- 
bilities resulting in the horizontal and normal branches 
only being present at the highest M), but the proposed 
scenario provides no explanation for the systematic sup- 
pression of horizontal and normal branches as 2-2.9 keV 
count rate drops. 

The above discussion has mostly focused on the role of 
M in the Z source tracks. As already mentioned in f|lj 
motion along atoll tracks appears to be mostly the result 
of changes in M . The spectral fits by LRH09 support 
this. In this respect, it is interesting to point to the 
atoll lower- banana branch in the HID of XTE J1701-462, 
which appears to be a continuation of the presumably M 
driven secular motion observed in the Z-like phase of the 
outburst. H owever, the presence of hysteresis in the state 
transitions ( Gladstone et al. | |2007) an d the phenomenon 
of parallel tracks ( [jVIendez et al. [11998 ) observed in other 
atoll NS-LMXBs suggest that a time avera ged response 
to ch anges in M might be involved as well (van der Klis 
20011). 



4.4. Kilohertz QPOs 

The luminosity range over which kHz QPOs are de- 
tected in XTE J1701-462 spans a factor of ~ 1 5-20, which 
is wider than seen in any other NS-LMXB (Ford et al. 
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20001. However, there is a fairly large gap in luminosity 
in wnich no kHz QPOs are detected, corresponding to 
selections G-K. In selections G and H the source shows 
Sco-like Z behavior, but is missing the horizontal and up 
per normal branches, where most of the kHz QP O detec 
tions in the Sco-like Z sources are made 



(see, e.g. 



van der 



Klis et al.|1997 Zhang etayi998| prnian et al':|20U2p . 
i'he behavior of A'i'E J 1701-462 in selections 1-K is sim- 
ilar to that of the bright atoll sources GX 9-1-9, GX 9+1, 
and GX 3-1-1, fo r which no kHz QPOs h ave been reported 
either (see, e.g.,|Wijna nds et al.|1998b| ). Surprisingly, we 
did not detect kHz QPOs m the Gyg-like Z selections, de- 
spite the fact that all thr ee persistent Cyg-like Z sources 



have shown kHz QPOs (Wijnands et al 



et al. 1998 



1998a Jonker 



2002). As pointed out in jJ3.1.1nowever 



selections A and B also show some differences from the 
Cyg-like Z sources in terms of CD/HID structure. Se- 
lection C is missing its horizontal branch, on which most 
kHz QPO detections were made in the Cyg-like Z sources. 
Finally, the fact that the Z-source-like and atoll-source- 
likc kHz QPOs in XTE J1701-462 are well separated in 
luminosity and have distinct properties (such as Q-value 
and rms amplitudes) suggests that they might be the re- 
sult of different excitation mechanisms, not only in XTE 
J1701-462 but also in other Z and atoll sources. 



4.5. Long-term modulations 
Quasi-periodic long- term variability has been observed 

tE 
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in nu merous LMXBs (|Charles et al 
2010 1 . The long-term modulations seen m the low-energy 
hght curve of XTE J1701-462 are closely related to the 
gradual changes in the CD/HID tracks of the source, and 
we therefore suggest that they are the result of changes 
in M . They are only seen at intermediate luminosities 
(not at the peak of the outburst or during the decay) . To 
study the time dependence of the long-term modulations 
we created a dynamical power spe ctrum of the long-term 
low-energy light curve (see, e.g., Clarkson et al."2003), 
by combining Lomb-Scargle periodograms ( [Scarglc 1982 ) 
of 75-day intervals, with the start dates of the intervals 
shifting by five days. The resulting dynamic power spec- 
trum reveals that the period of the long-term modula- 
tions varies between ^^20-50 days. Given the strong vari- 
ations in their period, the long-term variations are prob- 
ably not the result of orbital modulations. We also note 
that the sharp dips that are present near the maxima 
of the first few cycles are not related to the absorption 
events that are often seen in high-inclination LMXBs; 
here they are the result of spectral changes associated 
with the source entering the horizontal branch. 

The lack of a known orbital period for XTE J1701-462 
complicates the interpretation of the long-term modula- 
tions, although, given the brightness and duration of the 
outburst, the orbital period of XTE J1701-462 is likely 
to be on the order of days or longer (H07). Accretion 
disk precession as the result of radiation- induced warping 
has been suggeste d as a source for long-term modulation 



in some LMXBs ( Pringlel 1996[ [Wijers fc Pringle||1999, 
Ogilvie & Dubus 12001), although such modulations are 



usually discussed m terms of changes in viewing angle or 
obscuration (which we rule out as a so urce of subclass 
transitions in XTE J1701- 



disc could perhaps also result in actual modulations in 
M, as the result of varying torques from the secondary. 
Two NS LMXBs have shown modulations that may be 
relevant to the discussion of the long-term variations in 
XTE J1701-462. |Shih et al.| ( |2005| report on modula- 
tions with a period of ~40 days m the NS LMXB 4U 
1636-53. These modulations are accompanied by transi- 
tions between the island state and lower banana branch 



(see also Belloni et al.||2007 l and are therefore likely the 
result of modulations in M. Similar to XTE J1701-462, 
the modulations in 4U 1636-53 were not present at every 
luminosity level; they only appeared after the long-term 
average lum inosity of the source had shown a decrease. 
|Shih et al.| ( 2005 ) argue that such a decrease could lower 
the irradiation of the outer accretion disc, resulting in 
a drop in its temperature and possibly bringing it into 
a regime of unstable accretion that results in limit cy- 
cle behavior. Based on the length and luminosity of its 
outburst, XTE J1701-462 is expected to have an orbital 
period on the order of days or more (H07) , which is longer 
than 4U 1636-53 (3.8 hrs). Presumably, it therefore also 
has a larger accretion disc, which means that the onset of 
such an instability could already occur at a higher lumi- 
nosity, in line with what we observe for XTE J 1701-462. 
The long-term modulations in Cyg X-2 (Wijnands 
|et a l.||1996 Paul et al. 2000) are accomp anied by strong 
changes between various types of Z tracks ( Kuulkers et al. 
[T996J , similar to what we observe in XTE J 1701-462. 
Although the modulations in Cy g X-2 appear to be less 
regular than in XTE J1701-462, [Boyd fc Smale] ( |2004 1 
find that the lengths of the modulations cycles in Cyg 
X-2 are consistent with being integer multi ples (1-14) of 
the 9 .843 day orbital period of the system ( [Cowley et al " 
19791. 
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4.2 1. Although not 



explicitly discussed in these works, a precessing warped 



5. SUMMARY & CONCLUSIONS 

The outburst of XTE J1701-462 has revealed in great 
detail the evolution of a low-magnetic field NS-LMXB 
from near- or super-Eddington rates all the way down 
to quiescence. Most interestingly, XTE J1701-462 was 
observed to transform from a Cyg-like Z source to a Sco- 
like Z source at high luminosities, and from Sco-like Z 
to atoll source at low luminosties. A brief comparison 
with other NS-LMXBs suggests that the behavior ob- 
served in XTE J1701-462 is representative of the entire 
class of low-magnetic field NS-LMXBs. We conclude that 
these transformations between NS-LMXB subclasses are 
the result of changes in a mass accretion rate component 
M (which for most of the outburst is measured by the 
2-2.9 keV count rate) and argue that differences in neu- 
tron star spin/mass/magnetic field or viewing angle are 
of minor importance to the Z/atoU phenomenology. The 
evolution observed within the Z source phase of the out- 
burst shows a surprising evolution of the Z source tracks, 
with the horizontal, normal, and flaring branches disap- 
pearing one-by-one as M appears to decrease. This is 
at odds with scenarios in which the same M component 
increases or decreases monotonically from the horizontal 
to the flaring branch. We suggest that the M responsible 
for the subclasses does not vary significantly along the Z 
tracks, while it does along atoll tracks. Motion along Z 
tracks may instead be the result of instabilities, whose 
presence does depend on M. 
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